Detector resolution effects on hadronic mass moments in B — > Xlv 
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We show how to relate moments of the reconstructed hadronic invariant mass distribution to the 
moments of the physical spectrum in a model independent way. The only information needed on 
detector resolution functions are their first few moments. Theoretical predictions for the first three 
moments as functions of non-perturbative parameters are given for various lepton energy cuts. 
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Differential distribution for inclusive decay rates of 
heavy mesons can be calculated using the operator prod- 
uct expansion (OPE) [1]. The resulting decay distribu- 
tions are singular and can only be compared with ex- 
perimentally measured distributions after smearing them 
over a sufficiently large interval. For example, calculating 
the hadronic invariant mass distribution for the semilep- 
tonic decay B — > X c lv using the OPE, one finds to lead- 
ing order in A QCD /m b [2] 
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where s = m x is the square of the hadronic invariant 
mass, rnp = i(nT>D + 'imp*) is the spin averaged D-D* 
mass, and T denotes the total semileptonic width. Non- 
perturbative corrections to this differential rate yield 
even more singular terms. To compare this distribution 
with the smooth spectrum measured experimentally, one 
has to find observables which sufficiently smear the dif- 
ferential spectrum. The simplest such observables are 
moments of this spectrum, defined as 
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From Eq. (1) on obtains trivially 
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The corrections to the first two moments have been calcu- 
lated to order 1/to| and q. 2 /3q [2-4]. By measuring devia- 
tions from the above relation, one can therefore obtain in- 
formation about the size of the non-perturbative matrix 
elements. In turn, information on these non-perturbative 
parameters can be used in precise determination of CKM 
elements \V c b\ and \V u b\- 

The first two moments were first measured by the 
CLEO collaboration with a lower cut on the lepton en- 
ergy, to suppress the background in which the lepton 
originates from a semileptonic decay of the D meson [5] . 
More recently, the DELPHI collaboration presented mea- 
surements of these moments without a cut on the lepton 
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energy [6], while the BaBar collaboration measured the 
first moment as a function of the lepton energy cut [7] . 
This dependence on the lepton energy cut can be com- 
pared with predictions from theory, allowing for a more 
detailed test of the operator product expansion. 

In general, the reconstructed hadronic invariant mass 
spectrum dT/dsu is a convolution of the true spectrum 
dT/dsT and a detector resolution function P(sr—st, St), 
which describes the probability of measuring the recon- 
structed invariant mass sr, given the true mass St 
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where X denotes any charmed final state not equal to 
a single D or D*. In (4) we have allowed for different 
resolution functions for D, D* and X final states as is 
expected experimentally. Each term in this expression 
depends explicitly on the value of the lepton energy cut 
E cu t- In the remainder of this paper, this dependence 
will be suppressed, however all equations should be un- 
derstood as having this dependence. 

From the result of [2-4] for the first two moments of 
the differential decay spectrum 
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we can then obtain the spectra needed in (4) using 
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where r^oo = T D (,)(E cut ) are the decay rates for B — > 
D^lv in the presence of a lepton energy cut. Note that 
only for E QvA — can the values of T D (,) be obtained 
from the measured branching fractions. 

As explained above, the theoretical expression for the 
differential decay rate is given in terms of a singular ex- 
pansion. Thus, this spectrum can only be compared with 
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the experimentally measured spectrum after smearing it 
with some appropriate weight function. One can show 
that this weight function has to be a smooth function 
with width of order y/rrn, Aqcd- While the detector res- 
olution function does provide a smearing of the theo- 
retical spectrum, its width is only several hundred MeV 
[7]. Thus, additional smearing is required to compare the 
calculated spectrum with the measured spectrum. The 
simplest observables to compare are moments of the re- 
constructed hadron invariant mass spectrum. 

These moments can easily be obtained from the general 
expression given in Eq. (4) 
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Here we have defined 



P N {st) = Jds R (s R -m 2 D ) N P( SR -ST,s T ), 



(9) 



(10) 



where P stands for any of the three resolution functions. 
Thus, moments of the reconstructed spectrum are given 
by the true spectrum, weighted with a function which is 
determined by the detector response. 

The result in Eq. (6) can be further simplified by ex- 
panding the functions P N (st) in a Taylor series: 



P N (s T ) = J2(sT-mlrP, N 
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The coefficients P£f depend on the precise form of the 
detector resolution function, but we can obtain some es- 
timates using the expected form of the detector resolu- 
tion function. We expect that for a given value of st 
the function P(u), st) is peaked around u) = with some 
width a that may depend weakly on st- The resolution 
functions used by the Babar collaboration [7] have width 
of order a <~ 1 GeV 2 , while for a measurement using fully 
reconstructed B events we expect a ~ (700 MeV) 2 . We 
furthermore assume that the variations in the resolution 
function are slow compared to the width of the resolution 
function. If the width a varies by a small fraction e over 
the range of the hadronic mass, 
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then the coefficients P„ are suppressed by e for n > N: 
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We can now see how the expected form of the reso- 
lution function results in a simple organization of the 



expansion (11). Using the estimates (13) and ((st — 
™£>) n ) ~ Aqq D we estimate the different terms in 
Eq. (11) as follows: 
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Thus, if the resolution functions behave as expected, 
the terms with n > N + 1 in the relation (??) giving the 
measured moments in terms of the theory moments are 

A 2 

suppressed by e-^{ QC ° )- Note also that since present 

™b m b 

experiments have a substantially larger than Aqcd the 
contributions of low n theory moments to the measured 
iV-th moment are amplified by powers of a/ Aq CD . Ne- 
glecting the suppressed terms (n > N+ 1) we thus obtain 
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This is the main result of this paper and shows how mo- 
ments of the reconstructed hadronic invariant mass are 
related to the moments of the true mass distribution. 
The only ingredient in this relation are the first few mo- 
ments of the detector resolution functions. 

As an example of these simplifications we consider a 
simple model for the detector resolution function 
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Here, the width of the Gaussian distribution depends on 
the value of st, however the variation is slow compared 
to the width of the Gaussian itself 
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For this model the Taylor expansions of the first few func- 
tions P N , Eq. (11), are 
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The first two moments arc quite trivial, but the second 
displays already all of the features of the general organi- 
zation principles given in (13). 

To summarize, we give a detailed description of the 
steps necessary for a measurement of the first moment 
of the reconstructed mass distribution. The procedure 
for higher moments should be obvious from this. For 
each of the detector resolution functions, calculate its 
first moment for various values of the true invariant 
mass st- This gives the function P 1 (st)- Fit this func- 
tion to a third order polynomial in (s T — m 2 D ) and call 
the coefficients P\. The theoretical expression for the 
reconstructed moment is then obtained using Eq. (15) 
for N = 1 and using the theoretical expression for the 
((st — m 2 D ) n ) given in the Appendix. 

For many applications the detector resolution function 
P{sr — s t, st) does not depend on the value of the true 
invariant mass and is therefore only a function of the 
difference of the true and the reconstructed mass 

P(sr - st,st) = p(sr - st) ■ (19) 

In this case, the expression in Eq. (15) simplifies con- 
siderably, since the quantities P„ can now be expressed 
solely in terms of moments of the resolution function. We 
find' 

^=(^y-, (20) 

where p N = Jdss N p(s). This allows to write an even 
simpler relations between the moment of the recon- 
structed and true spectra. For example, for the first 



moment 

(s R -m 2 D ) = (s T -m 2 D ) + p\ (21) 
+ (Pd -Px) — + [Pd* ~ Px) — > 

Note that Eq. (21) is exact, i.e., it docs not contain un- 
known 6(tAq CD /toj corrections as Eq. (15) does. In gen- 
eral, for the N-th moment of the reconstructed spectrum 
the first iV moments of the resolution function, as well as 
the first N moments of the true spectrum are needed. For 
distributions with vanishing first moment, p 1 = 0, we see 
that the detector resolution function introduces no bias 
into the determination of the first moment of (st ~rn 2 D ). 
In practice, however, a detector always can miss particles 
and the first moments, p , are therefore expected to be 
negative. 

To summarize, we have shown how to relate moments 
of the reconstructed hadronic invariant mass spectrum 
to the moments of the true, calculable, spectrum, involv- 
ing only moments of the detector resolution functions. 
We have given expressions for the general case, in which 
the resolution functions can depend on the true invariant 
mass, as well as for the special case, in which the resolu- 
tion functions only depend on the difference between the 
true and the reconstructed invariant mass. 
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invariant mass moments. We furthermore would like 
to thank Oliver Buchmullcr, Hcnning Flachcr and Vera 
Luth for helpful discussion. This work was supported by 
the DOE under grant DOE-FG03-97ER40546. 



APPENDIX A: THE FIRST THREE HADRONIC MOMENTS 



In this appendix we give the expressions for the first three hadronic moments S n — ((st — ™f))™) [note that 
this definition of S n is slightly different than the one for A4 n in Eq. (2)]. Since the hadronic moments satisfy 
S n <~ (Aqco/ms)", the first three moments are non-vanishing to order l/m 3 B . The decay rate T, and the first two 
moments Si , S2 were calculated in the pole scheme as a function of the lepton energy cut in [2] . In [4] this calculation 
was repeated in several different mass schemes and explicit expressions were given for the dependence on the lepton 
energy cut. For completeness, we will repeat these results here and also give the result for the third moment S3, which 
has not been presented in the literature so far. All these results are given in the IS scheme, and for comparison with 
earlier works also in the pole scheme. They are written in the form 

S n (Eo) = S^(Eo) + S^(Eo) A + S^(E ) A 2 + S^(E ) A 3 + S^(E ) Xi + S^(E ) A 2 
(E ) Ax A + (E Q ) A 2 A + (E ) Pl + (E ) p 2 + (E ) 71 
+S( 12 ) (E ) T 2 + (E ) % + Si 14 ) (Eo) % + (E ) e + (E ) e| LM + S< 17 > (E ) eA , (Al) 

where E is the value of the lepton energy cut, and the coefficients S^(Eq) has the appropriate mass dimension. 
The parameter e = 1 and cblm = 1 denote the order in perturbation theory, for which we have used a s (m,b) = 0.22 
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TABLE I: Coefficients for Si(Eq) in the Poie scheme as a function of Eq. 
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TABLE II: Coefficients for S 2 (E ) in the Pole scheme as a function of E . 
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and Po = 25/3. The definition of the parameter A depends on the mass scheme used and is 



AP° lc = fn B mr' c + - Pl - T \- % (A2) 
2m b mi 
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The values for the coefficients S$ are shown in tables I to VI. The tables do not contain the eA term for non-zero 
lepton energy cut E . This is because the radiative corrections to the hadronic energy distribution with a lepton 
energy cut are not known; for details see Ref. [3]. In Appendix B we give analytical expressions required for the eA 
term with zero lepton energy cut. 



APPENDIX B: ANALYTIC EXPRESSIONS FOR PERTURB ATIVE CORRECTIONS TO NEW 

PARTONIC MOMENTS 

The leading pcrturbativc corrections to the hadronic invariant moments arise only from brcmsstrahlung contri- 
butions, since the virtual contributions are proportional to a (5-function which does not contribute to the moments. 
These bremsstrahlung contributions were calculated in [3] with a cut on the lepton energy. However, virtual diagrams 
do contribute to perturbative corrections to the terms supressed by powers of Aqcd/^6- Using 

sh = m% + 2m B m b (e - 1) + m\{\ + s - 2e ) . (Bl) 

the corrections proportional to a s AQCD/?nb can be obtained from the perturbative corrections to the partonic invariant 
mass so = Pc/ m b an d partonic energy eo = e c /mi,. For the first two hadronic invariant mass moments one needs the 
partonic moments M^°'°\ M^'°\ M^°\ M^ ' 1 ), M^ ^ and where we used the definition of [2] 

J ds de 

with m c = m c /mb. Perturbative corrections to these partonic moments can be obtained using the results of [8] and 
analytical results were given in [2] . Here we give the expressions for the additional moments needed to calculate the 
third hadronic invariant mass moments. Writing 
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TABLE III: Coefficients for S-i(Eo) in the Pole scheme as a function of Eq. 
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TABLE IV: Coefficients for Si(Eo) in the IS scheme as a function of Eq. 
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we find for the perturbative corrections 
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TABLE V: Coefficients for S^-Eo) in the IS scheme as a function of Eq. 





Si 


o<2) 


c.(3) 
a 2 


S 2 


S 2 


c( 6 ) 
S 2 


c( 7 ) 
a 2 


c.(S) 
a 2 


c( 9 ) 

b 2 


c (l0) 

6 2 


o(n) 

s 2 


s 2 


S 2 


s 2 


S 2 


S 2 


i, 2 





0.816 


3.188 


3.889 


1.73 


-1.95 


-1.28 


5.06 


-2.78 


-4.69 


0.67 


-0.06 


-2.38 


0.77 


0.1 


0.454 


-0.47 


0.356 


0.5 


0.795 


3.106 


3.784 


1.68 


-1.88 


-1.19 


5.08 


-2.55 


-4.75 


0.57 


-0.05 


-2.32 


0.76 


0.11 


0.469 


-0.465 




0.7 


0.77 


3.005 


3.651 


1.61 


-1.79 


-1.06 


5.13 


-2.23 


-4.83 


0.43 


-0.03 


-2.23 


0.76 


0.13 


0.499 


-0.458 




0.9 


0.738 


2.872 


3.471 


1.5 


-1.63 


-0.87 


5.24 


-1.74 


-4.99 


0.24 


-0.01 


-2.11 


0.75 


0.16 


0.562 


-0.446 




1.1 


0.702 


2.722 


3.258 


1.37 


-1.42 


-0.62 


5.48 


-1.08 


-5.24 


-0.02 


0.03 


-1.94 


0.74 


0.2 


0.678 


-0.43 




1.3 


0.668 


2.572 


3.031 


1.2 


-1.12 


-0.3 


5.94 


-0.21 


-5.66 


-0.34 


0.08 


-1.72 


0.74 


0.26 


0.888 


-0.411 




1.5 


0.641 


2.45 


2.826 


1.03 


-0.68 


0.11 


6.9 


0.98 


-6.41 


-0.7 


0.17 


-1.43 


0.76 


0.34 


1.302 


-0.388 




1.7 


0.638 


2.423 


2.746 


0.92 


0.13 


0.72 


9.27 


2.84 


-7.99 


-1.03 


0.35 


-0.94 


0.82 


0.47 


2.295 


-0.369 





TABLE VI: Coefficients for S3(Eo) in the IS scheme as a function of Eq. 



E Q 


<?(!) 
3 3 


S ™ 


5< 3 > 




4 5) 


<j(8) 

a a 




^ 3 


o(9) 


15 3 




„(12) 
S 3 


<j(13) 


e(14) 
a 3 


S ™ 


<j(18) 
15 3 


15 3 





0.87 


5.084 


11.102 


11.43 


-9.59 


-2.52 


-12.73 


-10.45 


10.06 


3.84 


-1.51 


-7.75 


1.08 


0.03 


3.094 


0.309 


3.441 


0.5 


0.84 


4.905 


10.7 


10.99 


-9.26 


-2.34 


-12.05 


-9.68 


9.54 


3.47 


-1.45 


-7.5 


1.06 


0.05 


3.067 


0.255 




0.7 


0.806 


4.7 


10.232 


10.47 


-8.85 


-2.12 


-11.14 


-8.7 


8.82 


3. 


-1.38 


-7.18 


1.04 


0.08 


3.1 


0.205 




0.9 


0.764 


4.445 


9.638 


9.78 


-8.29 


-1.81 


-9.8 


-7.31 


7.75 


2.35 


-1.28 


-6.76 


1.01 


0.12 


3.235 


0.148 




1.1 


0.719 


4.167 


8.977 


8.99 


-7.61 


-1.43 


-8.02 


-5.55 


6.3 


1.56 


-1.16 


-6.26 


0.98 


0.17 


3.558 


0.1 




1.3 


0.675 


3.896 


8.314 


8.15 


-6.83 


-0.96 


-5.68 


-3.38 


4.33 


0.62 


-1.02 


-5.7 


0.97 


0.25 


4.223 


0.062 




1.5 


0.641 


3.674 


7.739 


7.37 


-5.94 


-0.38 


-2.43 


-0.6 


1.52 


-0.47 


-0.84 


-5.09 


0.97 


0.35 


5.608 


0.028 




1.7 


0.631 


3.594 


7.486 


6.94 


-4.75 


0.44 


3.2 


3.45 


-3.39 


-1.76 


-0.58 


-4.35 


1.04 


0.52 


9.047 


-0.01 
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